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Abstract. We develop our method to prove quantum supcrintcgrability of an integrablc 
2D system, based on recurrence relations obeyed by the eigenfunctions of the system with 
respect to separable coordinates. We show that the method provides rigorous proofs of 
supcrintcgrability and explicit constructions of higher order generators for the symmetry 
algebra. We apply the method to 5 families of systems, each depending on a parame- 
ter k, including most notably the caged anisotropic oscillator, the Tremblay, Turbiner and 
Winternitz system and a deformed Kepler-Coulomb system, and we give proofs of quantum 
superintcgrability for all rational values of k, new for 4 of these systems. In addition, we show 
that the explicit information supplied by the special function recurrence relations allows us 
to prove, for the first time in 4 cases, that the symmetry algebra generated by our lowest 
order symmetries closes and to determine the associated structure equations of the algebras 
for each k. We have no proof that our generating symmetries are of lowest possible order, 
but we have no counterexamples, and we are confident we can can always find any missing 
generators from our raising and lowering operator recurrences. We also get for free, one 
variable models of the action of the symmetry algebra in terms of difference operators. We 
describe how the Stackel transform acts and show that it preserves the structure equations. 

Key words: supcrintcgrability; quadratic algebras; special functions 
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1 Introduction 

An ?i-dimensional quantum (maximally) superintegrable system is an integrable Hamiltonian 
system that not only possesses n mutually - commuting symmetry operators, but in addition, 
the Hamiltonian commutes with n — 1 other finite order partial differential operators such that 
the In — 1 operators are algebraically independent. We restrict to quantum systems of the form 
H = A n + V where A n is the Laplace-Beltrami operator on some real or complex Riemannian 
manifold and V is a potential function locally defined on the manifold. These systems, including 
the quantum anisotropic oscillator and the hydrogen atom have great historical importance, due 
to their remarkable properties, [U [2]. The order of a quantum superintegrable system is the 
maximum order of the quantum symmetries as differential operators (with the Hamiltonian 

*This paper is a contribution to the Special Issue "Symmetry, Separation, Super-integrability and Special 
Functions (S 4 )". The full collection is available at |http://www.emis.de/journals/SIGMA/S4.html| 
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excluded). However, we always choose the generators such that the maximum order is as small 
as possible. Systems of 2nd order have been well studied and there is now a structure and 
classification theory [31 HJ [5j [61 [8] . Until very recently, there were comparatively few known 
superintegrable systems of order > 3 and virtually no structure theory for the algebra generated 
by the symmetries. In the last three years, however, there has been a dramatic increase in 
discovery of new families of possible higher order superintegrable classical and quantum systems 
QUI EU HI Q3J Q31 [T51 [16] . The authors and collaborators, and others, have developed methods 
for verifying superintegrability of these proposed systems, [T71 [J3J HI EQl [21] . The first method 
developed by us, to verify superintegrability for 2-dimensional quantum systems, was based on 
a canonical form for symmetry operators of arbitrary order. This approach succeeded for several 
important systems, such as the caged anisotropic oscillator and the TTW potential for all ra- 
tional k = p/q, but it led to multi-term recurrence formulas for which very careful analysis was 
needed to verify finite dimensional solution spaces, [201 122j . Furthermore, the approach yielded 
no information about the structure of the algebras generated by the symmetries. In [22] the 
authors introduced a new method, based on the recurrence relations obeyed by the separated 
eigenfunctions of the Hamiltonian, and sketched its application to the caged anisotropic oscilla- 
tor. In this paper we further develop and make rigorous the special function recurrence relation 
method, and apply it to additional cases, including the TTW system. We show that this new 
method enables us to prove for the first time that the symmetry algebra for the TTW potential 
is closed for all rational k, and to compute the structure relations for the algebra. (In a series of 
recent papers [231 l24l l25l [26] Marquette has used a ladder operator method to determine higher 
order symmetry operators and structure equations for 2D superintegrable systems that separate 
in Cartesian coordinates. In particular, he has found the structure equations for the caged an- 
harmonic oscillator. Our method is related to his for separation in Cartesian coordinates, but 
more general in that it applies to systems that separate in any coordinate system, e.g., polar 
coordinates.) 

The basic idea that motivates the method is that, since the formal eigenspaces of the Hamilto- 
nian are invariant under action of any symmetry operator, the operator must induce recurrence 
relations for the basis of eigenfunctions corresponding to the associated coordinate system in 
which the eigenvalue equation separates. The recurrence relation method uses the known re- 
currence relations for hypergeometric functions in order to reverse this process and determine 
a symmetry operator from a suitable combination of recurrence relations. We can compute the 
symmetry operators and structure equations for the symmetry algebra by restricting ourselves 
to a formal "basis" of separated eigenfunctions. Then we appeal to our theory of canonical 
forms for symmetry operators to show that results obtained on restriction to a formal eigenbasis 
actually hold as true identities for purely differential operators defined independent of "basis" 
functions. (We have no proof that this reverse process will always work but conjecture that it 
will succeed whenever all the separated eigenfunctions are of hypergeometric type.) We start 
with a simple example on the 2-sphere, to introduce the theory. Then we consider an example 
in Minkowski space followed by a revisiting of the caged anisotropic oscillator. Finally we treat 
the TTW potential again. In four of these cases we give the first proofs of the closure and the 
structure of the symmetry algebras. 

2 A simple system 

Here we construct a proof of quantum integrability using recurrence formula techniques. As 
a trial system we consider a quantum Hamiltonian on the two-sphere: 



H = 8 2 e + cot 6d e + —TTpl + 



a 



(2.1) 



sin 



sin 



2 cos 2 k(p 
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where k = ~ and p, q are relatively prime positive integers. (Here, all our considerations are 
local. We do not require that the potential is globally defined on the two-sphere, or that there 
are any boundary conditions. Global issues can be examined on a case-by-case basis. Also, all 
parameters and variables can be complex, except for k which is rational.) This system is clearly 
integrable, since it admits a 2nd order symmetry, responsible for separation of the eigenvalue 
equation for H in spherical coordinates. Our aim is to show that this Hamiltonian admits an 
additional quantum symmetry, so is superintegrable. This is relatively straightforward to do 
using recurrence formula techniques. The typical separable solution in spherical coordinates has 
the form 



^ +|) (cos^)(cosV) 1/2 f/^(sin^), 



* = T, 

where T,U = P or Q are the solutions of Legendre's equation, a = k 2 (j — a 2 ) and ifj = kip. 
Indeed if we write = 0(6)&(ip) then these solutions satisfy the "eigenfunction" equations 
H^> = -n(n + 1)^ where H = L\ and L 2 ^> = -k 2 (N + ^) 2 ^, or 

/ k 2 (N + ±) 2 \ 
d 2 e + cot 6d e , 9 2 9 = -n(n + 1)0 



v sm 2 

and 

(«+^)— * h: 2 

where 

$ = e(0)$(<p) = rn (JV+ ^(cos0)(cosV) 1/2 ^(sin^). 

(Note that L 2 is the symmetry operator associated with variable separation in spherical coordi- 
nates.) Here no boundary conditions are implied and all parameters can be complex. To proceed 
further we observe the following recurrence formulas for Legendre functions 



D+(x)T?(x) = (1 - x 2 )^-T^(x) -{u + l)xTjf(x) = -{v - /i + l)T^ +1 (x), 
D-{x)T^(x) = (1 - x 2 )^(x) + vxT£(x) = (u + rfT^x) 



and 



C+(x)l?(z) = (1 - x 2 ) ll2 ^-T»{x) + _^ 2 l/2 Tt?(x) = -Ti+Hx), 

C-{x)T»(x) = (1 - x 2 f l2 ^{x) - _^ 2)X/2 mx) = {v + fi(v - m + l)rr X (x). 

These relations enable the shifting of the indices v and fi by ± 1. There are clearly similar 
relations for the functions Uj}{y). If we let N — > N + q then 

^ ^ r„ fc(7V+ ^ +P (cos0)(cosV) 1/2 ^ +g (sin^). 
We note that N can be quite arbitrary and now consider the function 

? + = T^ +p (x)U° N+q (y), 
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where x = cos 9 and y = simp. This function can be obtained from 

T = T^ N+ ^(x)U a N (y) 

via raising operators: 

A + (AOT = C+^^ix) ■ ■ ■ C^^D^y) • • • D»T 

= (-l)P+i(N-a + l) q T + . (2.2) 

Here (a)i = a(a+l) ■ ■ ■ (a+£— 1) is the Pochhammer symbol, and we assume that we have chosen 
bases for the separable solutions such that the same recurrence formulas and normalizations hold 
for all elements of the basis. Thus relation (|2.2p holds for all four elements in the basis. Similarly, 
we look at the possibility that N — > N — q. This can be obtained from To via lowering operators: 

A_(AT)T = C; (N+l) _ p+1 ( X ) ■ ■ ■ C- {N+i) {x)D- N _ q+l {y) ■ • • D N (y)T 

= (-l) p+q (-N-a) q (-n-k(N + ^ (n-k^N + ^j+lj T_, (2.3) 

where 

T- = T^ N+ ^ P (xm_ q (y). 

This result is independent of which basis we choose. 

Now consider the differential operator A = A + + A_. It follows from the form of the C 
and D operators that this is an even polynomial function of N + ^. This can be seen from 
the v —v — 1 symmetry of the operators D and the [i — > —fi symmetry of the operators C: 

C+(x) = CZ^x), D+(y) = DZ V _M- Note that (N + ±) 2 * = -k~\d% + ^g^)* for any 
eigenfunction $ of L2. It follows that under the transformation N — > — N — 1 we have A+ — > A_ 
and A_ — > A + . Thus, everywhere the term (N + ^) 2 ^ occurs in the expansion of A (where I 
is a positive integer) we can replace it by (k~ 2 L,2) and obtain a pure differential operator, 
independent of the parameters n, N. As a consequence we see that we have constructed a pure 
differential operator, which we also call A, and which preserves each eigenspace of H when 
acting on functions of the form 



T = Tn (N+1 *\x)U a N {y) = *(sin V0" 1/2 . 



Similarly, the operator A + — A_ goes to its negative under the mapping iV — > —N — 1, so it is an 
odd function of N+^. This implies that A = (l/(iV + i))(A + — A_) is an even function of N+^, 
so it can also be represented as a pure differential operator, independent of the parameters N, n 
and it preserves each eigenspace of H. We have constructed two partial differential operators A 
and A, each of which commutes with the Hamiltonian H on formal eigenspaces. Thus they 
act like symmetry operators. However, to prove this we must show that they commute with H 
when acting on any analytic functions, not just eigenfunctions. To establish this fact we use the 
canonical form for symmetry operators studied in our papers |20| [22] . 

Although we give the reasoning for this special system, our argument is quite general and 
immediately applicable to other systems. The operator L2 determines separable coordinates x, y 
for the system (in this case (x,ip). Now consider the commutator [H, A]. When acting on formal 
eigenfunctions To the commutator gives 0. We want to show that it vanishes identically. To do 
this we write [H, A] in canonical form by recursively replacing all second derivatives d% , dy in 
terms of H, L2 to obtain the expression 

[H, A] = A(x, y, H, L 2 )d 2 xy + B(x, y, H, L 2 )d x + C(x, y, H, L 2 )8 y + D(x, y, H, L 2 ). 
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This expression has to be interpreted as in [20], i.e., the parameters H, L 2 must be moved to 
the right before being identified as operators. Now applying this operator to any eigenfunction 
we obtain 



[H, A] T = (Ad 2 xy + B8 X + Cd y + D)T = 0, 



(2.4) 



for all choices of the parameters H, L 2 . Noting that we have 4 linearly independent choices for 
To = Pi{ x )Qj{y)i we can write (|2.4p as a set of 4 homogeneous equations for the 4 unknowns 
A, B, C, D: 

( P[(x)Q[(y) PiWQtiy) Px{x)Q' 1 {y) t\{x)Q x {y) 
P{(x)Q' 2 (y) P[(x)Q 2 (y) Pi{x)Q' 2 {y) P 1 (x)Q 2 (y) 
P!>{x)Q'M P^x)Q 1 (y) P 2 (x)Q' 1 (y) P2(x)Q 1 (y) 
\ P^(x)Q' 2 (y) P^(x)Q 2 (y) P 2 (x)Q' 2 (y) P 2 (x)Q 2 (y) 

It is convenient to introduce the determinant function 



\ 


< A \ 








B 









C 







/ 


\ D J 




W 



W(h(x), gi (y),f 2 (x),g 2 (y)) 



fl(x)g' 2 (y) 

fi(x) f 2 (x) 
fi(x) &(x) 



f[(x) gi (y) 
fl(x)g 2 (y) 

f 2 (x)gi(y) 
f' 2 {x)g 2 {y) 

9i(y) 
g'i(y) 



h(x)g[(y) 
fi(x)g' 2 (y) 
h(x)g' x {y) 
h{x)g' 2 {y) 

92(g) 

92(y) 



h{x)gi(y) 
fx(x)g 2 (y) 
h{x)g\{y) 
f2(x)g 2 (y) 



We then note that for our particular system 
W(P 1 ( a: ),P 2 ( a; ) ) Q 1 (y),Q 2 (y))^0, 



except at isolated points, since the Wronskian of two independent solutions of a separated 
eigenfunction equation is nonzero. Thus we conclude that A = B = C = D = 0. Consequently 
[H, A] = identically. This proves that A = L3 is a symmetry operator for the system. The 
same proof shows that A = L4 is also a symmetry operator. Since both L3 and L4 fail to 
commute with L 2 , each must be algebraically independent of H = L\,L 2 . Thus this system is 
superintegrable for all rational k. 

Example 1. If k = p = q = lwe have the familiar superintegrable system on the 2-sphere 
where the corresponding Schrodinger operator is 



H = di + cot 9dg + 



1 



sin 



-d 2 + 



o 



sin 2 9 cos 2 ip ' 



i.e., the Laplacian on the two sphere s 2 + s\ + s 2 



potential V(s) 
and L4 



r . This corresponds to system 



3 — 1 (in Cartesian coordinates Sj) plus the 
in the list @]. Then we have L3 = A + +A_ 



(^I) (A - 



A_) with A+ and A_ where 



A + (N)f = ((1 - y 2 )d y - (N + l)y) ( VT^d x + {N + * )x 



f 



x" 



(cos^fy - (A + 1) sin-i/O -d e + 



A_(A)/ = ((1 - y 2 )d y + Ny) (^/T 



A + - ) cot 1 



x 



2 a ; 



{N + ±)x 

VT^ 2 



f 
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= (cos ^ + N sin i/j) {-d e - [n + ^ cot e) f, 

where x = cos 9 and y = simp. Note that in this example tp = <p since k = 1. 
We choose the standard spherical coordinates on the sphere viz. 

si = sin 9 cos ip, s 2 = sin# simp, S3 = cos 9 
and the corresponding expressions 

Ji = s 2 d S3 - s 3 d S2 , J 2 = s 3 d Sl - sid sa , J 3 = si<9 S2 - s 2 d Sl . 
The Hamiltonian can then be written as 



H = J( + Jf + Jf + ^. 
From these calculations we deduce that 

S2£3 



o 



L A = -2J 1 , L 3 = {J 2 ,J 3 } + J 1 + 2a- 



where the separation equation in polar coordinates is due to L 2 = J 2 + a/cos 2 tp with eigenvalue 

2' 



-(N + i) 2 . (Here, {A,B} = AB + 5^.) We also note the relation 2L 3 + L 4 = [L 4 ,L 2 ]. Thus 



we see that [L4, L 2 ] is not proportional to L 3 . 

We can compute the general structure relations for the symmetries of system (]2.ip . We obtain 
[L2, L4] = -2qk 2 L 3 - q 2 k 2 L A , 

[L 2 , L 3 ] = -q 2 k 2 L 3 + 2qL 4 L 2 = q 2 k 2 L 3 + q 3 k 2 L 4 + q{L 2 , L 4 }. 
To compute [Z/3,L4] we need some preliminary results. We make note of the identities 
A + (N - e?)A_(iV)T = (-l)«(a - N) q (-N - a) q 

-n-k{N + ^\\ L-k\N+^\+l \ T = F 1 (n,N)To, 
A_(iV + q)A + (N)T = (-l) q (N-a + l) q (N + a + l) q 

-n + k[N + ^J (n + k^N+^+l) T = F 2 (n, N)? . 

Recall that E = — n(n+l) is the eigenvalue of the Hamiltonian and — k 2 (N+^) 2 is the eigenvalue 
of L 2 corresponding to a basis function. Using the property (— a)i = (— 1)^(q — I + l)i we note 
that Fi(n,-N - 1) = F 2 {n,N) and Fj(-n -1,N) = Fj(n,N) for j = 1,2. Consequently, 
F + = Fi(N,n) + F 2 (N,n) is an even polynomial function in N + | and a polynomial function 
of n(n + 1), hence when acting on separated basis functions F + = P + {H,L 2 ) is a polynomial 
in H and L 2 . In fact, F + = P(H, L 2 ) as pure differential operators, independent of basis. The 
proof of this fact is analogous to that given above. We write the operator F + — P + in canonical 
form Ad 2 y + bd x + Cd y + D . Then on an arbitrary eigenbasis we have 

(Ad 2 y + B8 X + Cd y + D)^ n = 0. 

It follows via the usual Wronskian argument that A = B = C = D = 0. Thus F + — P + = 
identically. Similarly F~ = — ^ {F x (n, N) -F 2 {n,N)) is an even polynomial function in N + t> 

N + j 

and a polynomial function of n{n +1), so F~ = P~(H, L 2 ) is a polynomial in H and L2. 
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[L 3 ,L 4 ] =qLl-2P-(H,L 2 



Now it is straightforward to obtain 

There is of course an extra constraint. In fact 

L 2 L 2 = -k 2 L 2 + qk 2 L±L 3 + 2k 2 P+(H, L 2 ), 
and, symmetrizing, we find 

{L 4 , L 4 , L 2 } = -6k 2 Ll - q 2 k 2 Ll - 3qk 2 {L 3 , L 4 } - Wqk 2 P~(H, L 2 ) + 12k 2 P + (H, L 2 ), 

where {A, B,C} is the symmetrizer of 3 operators. The structure is clearer if we define R = 
-2qk 2 L 3 - q 2 k 2 L 4 . Then we can rewrite the structure equations as 

[L 2 , L 4 ] = R, 

[L 2 ,R] = -2q 2 k 2 {L 2 ,L 4 } - q A k A L 4 , 
[L 4 , R] = 2q 2 k 2 L\ - 4qk 2 P~(H, L 2 ), 

+ {L 4 , L 4 , L 2 } - 1—L 2 - 12k 2 P + (H, L 2 ) + Wqk 2 P~(H, L 2 ) = 0. 
2q z k z 2 

This shows that the symmetry algebra is generated by the symmetries H, L 2 , L 4 and is closed 
under commutation. 

Example 2. We consider our system for the case k = ^, i.e., p = 1, q = 2: 
V 



2^Js{ + s'i(si + y/s'{ + s'i) 4sin 2 <9cos 2 (^)' 

where we now use the fact that ip = ip/2. This is a special case of system [S7] in [3]. We form 
the functions A + and A_ as before: 

A + (N)f = (cos yd,/, - ( A + 1) sin ip) (cos ipd^ - (A + 2) sin-i/0 



A_(A)/ = (cosy^ + AsinV)(cosya v , + (A r -l)sinV) (-fy - ^ ^A + cot^ /. 



The Hamiltonian becomes 

H = d 2 + cot ^ + — ^- 5^ + " . 

4 sm v 4 sm p cos z tp 

We obtain 



53(24 + 52-2^^ + 5!) 
2V s i + s : 



L 4 = -2{Ji, J 3 } + a cot0 tan 2 V = -2{Ji, J 3 } + a 

s,vsf + s 2 



L3 = — { J3, J3, J2} — L 4 + 2 cot ( a tan y H — sin ip cos ^ ) 



,2 .,. , o„„„o„/.\ \ fl , 4- a ( 2 



+ 2 atan 2 V + -(1 - 2cos2V>) <9 e + a cot 



3 J \ cos 2 tp 

The operator describing the separation in polar coordinates is 

L 2 = Ud 2 + 



^ cos 2 ip t 

which corresponds to eigenvalue —\{N + \) 2 ■ The basic commutation relation is L3 + L 4 
[U,L 2 ]. 
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3 Another system 

We can further investigate these ideas for the case of the special potential 

(x + iyf 

V = a-^ ^-r, (6.1) 

(x 2 + y 2 ) 4 

real in Minkowski space with coordinates x\, X2 where x = xi, y = 1x2 ■ In polar coordinates 
the Schrodinger equation has the form 

d 2 + -d r + \d 2 e + ^e m - e] tt(r, 9) = 0. 



If we write E = — (3 2 then typical separable solutions in polar coordinates are 

Cn(/3r)C n/3 {5e 3ie ), 

where C u (z) is a solution of Bessel's equation, = R(r)@(9) and a = —95 2 . Indeed 

(d 2 e - 95 2 e eie + n 2 )e(8) = o. 

We now construct operators which induce the transformation Q — > £1 ± 3 on the basis func- 
tions ^q. To do this transparently we use the variable w = 5e 3ie . We make the observation 
that 

$+ _ (_ 8r + "±») (_ a . + fl±i) (-a. + 5) (-a„ + : ^* n+3 . 

Similarly we note that 

...( ft + £zi)(* + 2zl)(* + 9) * n ^ n _ 3 . 

Clearly $ = <£_)_ + is an even function of 0, hence interpretable as a pure differential 
operator. Similarly the operator <3?_|_ — <£_ is an odd function of f2, so = — is 

a pure differential operator. We deduce as previously that [H, <£] = [iJ, 5>] = 0, hence we have 
constructed a quantum superintegrable system. 

We can extend these ideas to consider the potential 

= a (x + iy) k ~ 1 
(x — iy) k+1 ' 

where k = p/q and a = —k 2 5 2 . The solutions have the form 

^n = Cn(Pr)C qn/p {5e^^). 

In order to map solutions of a fixed (3 eigenspace into solutions we can use the transformations 
Q — > £1 ± p. These transformations can be performed by the differential operators 



and 



r J \ pw J 



-8 W + ^ + q - 1 ) : ->• /3 P ^ +P , 



.__ (aH .°±iz£y..k + fT) f 8 „ + ^ 



r J \ pw 
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x ( d w + ^ + 1 - q) : *n -> ^n-p, 
\ pu> / 

where w = 5 exp(ip9/q). If we make the transformation 0, — > —Q then we see that <£_(— f2) = 

(_i)p+«$ + (n). 

Consequently there are two cases to consider: 

(a) p + g even: $(+) = $ + + $_ is an even function of f2, hence is a pure differential operator, 
and <J>(~) = — 3>-) is also a pure differential operator. 

(b) p + (7 odd: <3?( + ) = $ + — <J>_ is an even function of f2, hence is a pure differential operator, 
and <J?(~) = ^($+ + $-) is also a pure differential operator. 

Thus we have superintegrability in both cases. 

Further, we can prove the finite closure of the symmetry algebra and construct the structure 
equations. We write 

L 3 = $ (+) , L 4 = $H 

and let L 2 be the differential operator whose eigenvalue corresponds to Q 2 , i.e. 

L 2 e(9) = {-8 2 + ^S 2 e 6 A 6((9) = n 2 Q(9). 

Then a direct computation verifies that the symmetry algebra structure relations are 

[L 2 , L 4 ] = R, [L 2 , R] = 2p 2 {L 2 , L 4 } - p A U, [L 4 , R] = -2p 2 L 2 , 
together with the Casimir 

R 2 - ^{L 2 ,L 4 , L 4 } + jp A L 2 + l6p 2 H 2p = 0. 

Here, {A, B} = AB + BA, {A, B,C} is the analogous 6-term symmetrizer of 3 operators, and 
R = 2pL 3 + p 2 L 4 . These relations hold no matter whether p + q is even or odd. 

This shows that the symmetries H , L 2 , L 4 generate the symmetry algebra, and that it closes. 



Example 3. Take k = p = q = 1. Then 



<5 2 <5 2 e 2i9 w 2 



(x — iy) 2 r 2 r 2 ' 

and 

*+/ = (-dr + -1 (-a* + -) f, *_/ = f 9 r + - + - ) /• 



r / \ to / \ r J \ w 

The Hamiltonian is 

H = d 2 + a 2 



y (x — iy) 2 

This corresponds to system [E14] in [3]. We find 

2 i 1 5 2 w 

Li = ~^{d x ~ id y ), L 3 = --{xd y - yd x ,d x - id y } + -{d x - idy) + 2— 

and the constant describing the separation of variables in polar coordinates is 

L 2 = w 2 d 2 v + wd w + w 2 = -{xdy - yd x ) 2 + 5 2 - + %V , 

x — iy 

which corresponds to the eigenvalue Q 2 . We also have the structure relation [L 2 , L 4 ] = 2L% + L 4 . 
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Example 4. Take k = p = 2, q = 1 and w = 5e 2id . Then 
(x — iy) 3 r 2 

and 



2w 

I "I ' 1 "HI ' ri 

r J \ r J \ Aw J 



The Hamiltonian is 



]2 , a2 ax2 x + l V 



H = di + di - 4<5 2 



x ' v {x-iy) 3 ' 
This corresponds to system [E8] in We find 



L 3 = -— {xdy - yd x ,d x - id y , d x - id y } + L 4 - 4-^ (2(rd r + wcta) + 1, 

1 .„,, .11! 



L 4 = -=(9,-^ + 4-2, 
and the constant describing the separation of variables in polar coordinates is 

L 2 = 4(w 2 dl + wd w + w 2 ) = -{xd y - yd x f + 45 2 (Z±M) * , 

which corresponds to the eigenvalue O 2 . We also have the structure relation [L2, L4] = ^(L^+Li). 

4 The caged anisotropic oscillator revisited 

In [22] we introduced the recurrence relation method by sketching a proof that the caged 
anisotropic oscillator was quantum superintegrable. Here we will provide more details and 
show that the symmetry algebra always closes. This result is not new [23] but we include it here 
to illustrate explicitly how it is obtainable from recurrences obeyed by Laguerre functions. The 
system is 

I_ n 2 I_„2 

it o2 1 q2 2 2 2 2,4 "1,4 "2 

H = d x +d -fi ± x -i4y + + o ' 

where ^\ = pfi and ^2 = and p, (? are positive integers that we assume are relatively prime. 

We look for eigenfunctions for the equation H^> = X^/ of the form ^ = XY and find the 
normalized solutions 

Xn = e-^x^+hL^^x 2 ), Y m = e-^y^L^^y 2 ), 

where the L"(x) are associated Laguerre functions [27]. Separation in Cartesian coordinates is 
determined by either of the symmetry operators 

T fl2 2 2,4 "! t fl2 2 2,4 2 

x y 
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where H = L\ + L 2 . For the separated solutions given above we have the eigenvalue equations 
Li^! = \ x fy, L 2 fy = X y ^, where 

X x = -2\x\{2n + a\ + 1), X y = -2p 2 (2m + a 2 + 1). 

Thus Hty = E^> where the energy eigenvalue is 

E = -2p(pn + qm+ pa\ + p + qa 2 + q). 

In the foregoing we impose no boundary conditions and the Laguerre functions are stand-ins 
for either of the two linearly independent solutions of the second order ordinary eigenvalue 
equations. In particular, n, m are allowed to be complex. In order that the eigenspace of H 
with eigenvalue E remain invariant under the action of a recurrence operator that changes m 
and n we must keep pn + qm constant. One possibility that suggests itself is that n — > n + q, 
m — )• m — p. A second possibility is n — > n — q, m — >■ m + p. 

Now note the recurrence formulas for Laguerre functions (or confluent hyper geometric func- 
tions) viz 

z^L%{z) = pL a p {z) - (p + a)^-i(^) = (p + 1)^+1(2=) - (p + 1 + a - z)L a p {z). 

We apply these for the cases z = pix 2 , z = u 2 y 2 and make use of the eigenvalue equations 
for X x , X y , respectively. Again, we can choose two linearly independent solutions for each 
eigenvalue equation, each of which satisfies this recurrence. Then, considering the symmetry 
operators as acting on basis functions ty n = X n Y n , we have the recurrences 

D + (m,x)X n = ^ - 2xp A - m + p^ x 2 + X n = -A t i 1 (n + l)X n+1 , (4.1) 

D~(m,x)X n = (^d 2 + 22^10* + Pi + Ax 2 + 4 x2 Gl ^j X n = -4pi(n + ai)X n _i, (4.2) 

D+(u 2 ,y)Y m = [d 2 y - 2 W2 d y - p 2 + p|y 2 + ^p^j Y ™ = -4p 2 (m + l)Y m+1 , (4.3) 

D-(fi 2 ,y)Y m = (^d 2 + 2 Vi i 2 d y + p 2 + ply 2 + ^p^j Y m = -4p 2 (m + a 2 )Y m _ 1 (4.4) 

for either basis solution. Writing u = m + kn where k = p/q, we have m = u — kn and we can 
characterize a formal eigenfunction corresponding to energy E = — 2u(qu + pa\ + p + qa 2 + q) 
by * n = X n Y m . Let = D+(pi, x)«D-(p, 2 , yf and $~ = D~(u 1 ,x) ( 'D+{fi 2 ,y)P. By direct 
calculation, using recurrences (I4.ip - (|4.4p . we verify the relations 

= (-4p 1 ) 9 (4p 2 ) p (n + l) q {-u + kn- a 2 ) p ^ n+q , (4.5) 
= (4pi)«(-4p 2 f (-n - ai ) q (u - kn + l) p ^ n - q . (4.6) 

(Note that this system is very simple to analyze compared to the other systems we study in 
this paper because these two operators are defined independent of n. Hence by the arguments 
that we have given for previous examples, each is a symmetry for the caged oscillator system 
that is algebraically independent of the pair L\, L 2 . No symmetrization or antisymmetrization 
is needed.) Thus the system is superintegrable. Now we construct the operators 

®i(n) = <3?+(n - g)$~(n) : 
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* n — ► (-16M?) 9 (-16/il) p (n -q + l) q (-« + fe(n - 9) - a 2 ) P (-n - - fen + l) p * n , 

<J> 2 (n) = 3>~(n + g)<E> + (n) : 

* n — > (-16M?) 9 (-16/il) p (-n - g - a x ) q (u - fe(n + q) + l) p (n + l) g (-« + fen - a 2 ) p * n - 

Though we have indicated a dependence of operators $1, $2 on n in order to compute their action 
on a formal eigenbasis, in fact we see from relations (|4.ip — (|4.4jl that they are pure differential 
operators, independent of the parameter n. Further they commute with both L\ and L 2 . Hence 
by an argument that we have given for a previous example, they must be polynomials in the 
symmetries L\ and H. Making the replacements 

H + 2fi(pai +p + qa 2 + 9) L\ — 1u,\{a\ + 1) 

u «->■ , n f-> 



2fiq A\x\ 
we find 

* a = (-16m?) Vie^ ( Ll " 2 ^ 1 + 1) - <Z + 1) 



if + 2/i(pai + p + ga 2 + q) ( L\ - 2\i\ {a\ + 1) 

+ fe ; 9 — «2 



2/ig V 4 A*i 

Li -2 m (ai + 1) 



4/x 



ai 



g + 2/i(pai + p + ga 2 + 9) L x - 2^i (01 + 1) 
2^9" 4^ + 1 

$ 2 = (-16a?)'(-16uir (~ Ll " 2 ^ 1 + 1) - 9 ~ ai) 

# + 2/x(pai +p + ga 2 + q) , ( L\ - 2//i(ai + 1) 

fe ; 1-9+1 



2"9 V 4/ii / p 
Li- 2^i(oi + l) | 1 

/ q 

H + 2/i(pai + p + 902 + 9) , , Li ~ 2/xi(ai + 1) 

+ fe "j «2 W„ = P 2 (ii, Li)V n . 



2/i9 4//i 



The constant terms in the expansions of these relations should be interpreted as the constants 
times the identity operator. Using these results and the eigenvalue formulas L\^ n = —2p\{2n + 
a\ + l)^ n , L 2 ^ n = — 2fi 2 (2u — 2fen + a 2 + l)\l/ n we can derive the structure equations for the 
symmetry algebra by acting on a formal eigenbasis. Then we can use our previous argument 
to show that the structure equations must hold independent of basis. Let L 3 = <£ + + <$~ and 
L 4 = <1> + — Then we have 

[Li, L 3 ] = -Ap,pqL A , [Li,L A ] = -4ppqL 3 , 

[L 3 ,L 4 ] = -2P 1 (H,L 1 ) + 2P 2 (H,L 1 ), L\ = Ll + 2P\(H, Li) + 2P 2 (H, L\). 

Thus the symmetry algebra closes. We can take the symmetries H, L\, L 3 as the generators 
with R = [L\, L 3 ] = —AppqLi and rewrite the structure equations as 

[L 1 ,L 3 ]=R, [L 1 ,R] = 16p 2 p 2 q 2 L 3 , [L 3 , R] = 8fipqPi(H, L x ) - 8fipqP 2 (H, L x \ 

*p2 _ T 2 



I6p 2 p 2 q 2 



R = L% — 2P l (H, L x ) - 2P 2 (H,L 1 ). 
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Example 5. We take the case of equal frequencies: H\ = fj, 2 = f^, so p = q = 1. The correspon- 
ding recurrence operators are 

® + f = (dl - 2 Xf id x - ix + ii 2 x 2 + ^\ (d 2 + 2yiid y + n + i?y 2 + ^\ f, 

$-/ = f aj - 2y/i^ -n + iJL 2 y 2 + J f<9 2 + 2x/i9 x + /x + ,uV + /, 

where Aj = j — a 2 . To proceed we need to have available the following operators: 



L x = d 2 x - fx 2 + 4, L 2 = d 2 - ii 2 y 2 + M = (xd y - yd x f + A^ + A 2 ^. 



where Li, L 2 come from our general theory and we recall that H = L\ + L 2 . We calculate the 
symmetry operators that our method implies. 

L 3 = $+ + <£>_ = 2L X L 2 + 4/i 2 M - 2/j 2 , L 4 = = n[M, L ± ]. 

Then we have 

[Li, L 3 ] = -4//L 4 , [Li, L 4 ] = -4/iL 3 . 
This method implies the existence of the symmetry operator M from the expression for L3. 

5 The TTW system 

A similar but more complicated procedure works for the quantum TTW system |12( [13] . Here 
the Hamiltonian is 



H = d$ + -d r - wV + —[81+ , + 



r 2 V sin 2 (£;0) cos 2 (^) y ' 
where we take = | as before. The general solution of the eigenvalue problem H^f = E^f is 

$r = e -f r V( 2n+a+6+1 )L^ 2n+a+b+1 )(^ (5.1) 

where we have taken a = k 2 (^ — a 2 ) and j3 = k 2 (j — b 2 ). The L-functions are associated Laguerre 
and the P- functions are Jacobi, not polynomials in general, [27]. We consider the functions 

n = e -7 r \^ 2n+a+b+ ^L k ^ 2n+a+b+1 \ur 2 )P^ b (x) = Y^(r)X^ b (x), 

where x = cos(2k0), X = P or Q and Y = S or T. (Here n is obtained from $ by a gauge 
transformation to remove the angular factors (sin(fc#)) a+ 2 (cos(k6)) b+ 2 .) By this we mean that 
Pn' b {x) is a Jacobi polynomial if n is an integer. Otherwise it is given by its hyper geometric 
expression. If X = Q then this denotes the associated second solution of the Jacobi differential 
equation. Similar remarks apply to the choice of Y = S. We have defined this function to be 

S£(r) = e -^\k{2n + a + b + l) L k[2n + a + b+l)^ 

and T^{r) to be a second independent solution. The energy eigenvalue is given by 

E = -2oj (2(m + nk) + 1 + (a + b + l)k) (5.2) 
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and A = k{2n + a + b + 1). The separation equation for Q(6) is 

L * B = { 9 > + *m + dk) m = - fc2(2n + ° + 6 + 1)20(9) = - A2Bm 

and L2 is a symmetry operator for the system. Under the gauge transformation L2 goes to 
a symmetry that we shall call L2 and which has the same eigenvalues. We see from the expression 
for E that in order that an energy eigenvalue be unchanged for different values of m, n we must 
fix u = m + nk. The two transformations 

n — > 11 + q, m — >• m — p and n — > n — q, m — >• m + p 

will each achieve this. 

Consider the functions Xn' h (x). If we want to raise or lower the index n we can do so with 
the operators [27] 

J+X% b (x) = (2n + a + 6 + 2)(1 - x 2 )<9 x X*> fe (x) 

+ (n + a + 6 + l)(-(2n + a + 6 + 2)x - (a - b))X% b (x) 
= 2(n + l)(n + a + b+ l)X^i{x) 

and 

J~X^ b (x) = -(2n + a + 6)(1 - x 2 )d x X%\x) 

- n((2n + a + 6)2; - (a - = 2(n + a)(n + 6)X^ 1 (x). 

Similarly, for the functions y^(R) = co A ^ 2 Y A (r) where R = r 2 we can deduce the relations [27] 

KX m yi(R) = \{A + 1)8 R - 1 - ^(,4 + 1)} y£(R) = -u>y£+l(R), 

Kl, m y£W = {(~A + l)d R - | + ±A(1 - A)\y*(R) = -cu(m+l)(m+A)y£- 2 (R). 

(Note that E = -2w[2(m + nk) + 1 + (a + b + l)fc], A = A;(2n + a + b + 1) and, effectively, the 
operator K + is lowering m by 1 and raising n by g/p, whereas the operator K~ is raising m 
by 1 and lowering n by q/p. Here, is fixed.) We now construct the two operators 

"+ = K A+2{p-l),m-{p-l) ' ' ' K A,m J n+q-l " ' J n ( 5 - 3 ) 

and 

= ^A-2(p-l),m+p-l ' ' ' -^A.m^n-g+l " ■ 

When applied to a basis function ^ n = 3^m(^)^« f° r fixed u = m + fcn, (so m = u — kn) 
these operators raise and lower indices according to 

S+* n = 2"(-l) p o; p (n + l) q (n + a + b+ l) q V n+q , (5.5) 
S_^ n = 2 q u p (-n - o) ff (-n - b) q (u -kn + l) p (-u - k(n + a + b + l)) p 4v_ 9 , (5.6) 

where (q) 9 = (a) (a + 1) • • • (a + 5 — 1) for nonnegative integer q, and we note the relation 
(— a)g = (— l) q (a — q + l) q . From the explicit expressions (|5.3[) . (j5.4|) for these operators it is easy 
to verify that under the transformation n — > —n—a—b—1 we have — > H_ and H_ — > H + . Thus 
E = H + + E_ as a polynomial in n and it is unchanged under this transformation. Therefore it is 



A Recurrence Relation Approach to Higher Order Quantum Superintegrability 



15 



a polynomial in (2re + a + 6+l) 2 and u. As a consequence of the relation A = — k 2 (2n + a+b+l) 2 , 
in the expansion of S in terms of powers of (2n + a + 6+l) 2 and E we can replace (2n + a + 6+l) 2 
by L2/k 2 and E by H everywhere they occur, and express H as a pure differential operator, 
independent of parameters. (Note that in the expansion of S in terms of the parameters, a 
term WE is replaced by WH with the W operator on the left.) Clearly this operator, which we 
will also call E commutes with H on the eigenspaces of H. However, the same argument as used 
in Example Q] shows that in fact S commutes with H in general, thus it is a symmetry operator 
for H. 

We can also easily see that under the transformation n — > — n — a — b — 1 the operator 
H + — H_ changes sign, hence the operator 5 = (l/(2n + a + b + 1))(H + — 5_) is unchanged 
under this transformation. Again, making the replacements (2n + a + b + l) 2 by L2/A; 2 and tt 
by —2 — (iJ + 2uik{a + 6 + l))/4u; we can express H as a pure differential operator, independent 
of parameters, and it is a symmetry operator for H. Each of these symmetries has a nonzero 
commutator with L2, so each is algebraically independent of the set H, L<i- This proves that 
the TTW system is quantum superintegrable for all rational k. We set L3 = S, L4 = S. 

Using the explicit relations (I5.5p . (I5.6f) for the action of the raising and lowering operators S± 
on a basis we can obtain very detailed information about the structure of the symmetry algebra 
generated by L2, L3, L4. Applying the raising operator to a basis function, followed by the 
lowering operator, we obtain the result 

~_(n + g)3+(n)¥ n = {-l) p A q uj 2p {n + l) q (n + a + l) q (n + b+ l),(n + a + b + l) q 

x (-u + fcn) p (it + k(n + a + b + 1) + l) p ^ n = ^ n ^ n . (5.7) 

Reversing the order we find 

S+(n - q)E_(n)^ n = (-l) p 4 q co 2p (-n) q (-n - a) q (-n - b) q (-n -a- b) q (u -kn+ l) p 

x (-u-k(n + a + b+l)) p ^ n = r] n ty n . (5.8) 

Thus the action of the operator H_(n + g , )5+(n) + H + (n — g)H_(n) on any basis function fy n 
is to multiply it by £ n + rj n . However, it is easy to check from expressions ([5.5p . fj5.6[> and 
from (I5.8p . (|5.7p that under the transformation n — > —n — a — b — 1 we have 3_(n + g)H + (n) <-> 
H + (n — g)S_(n) and ^ n -H- ry n . Thus = H_(n + g)3 + (n) + H + (n — g)H_(n) is an even 
polynomial operator in (2n + a + b + 1), polynomial in tt, and ^ n + % is an even polynomial 
function in (2n + a + 6 + 1), polynomial in u. Furthermore, each of H_(n + g)H + (n) and 
H + (n — g)H_(n) is unchanged under the transformation u — > —u — (a + b + 1) — 1, hence 
each is a polynomial of order p in {2u + (a + b + l)k + l) 2 = E 2 /Aco 2 . Due to the multiplicative 
factor ui 2p in each of these expressions we conclude that is a symmetry operator whose 
action on a basis is given by a polynomial operator P^ + \H 2 , L2, w 2 , a, 6). In fact, 

H(+)=p(+)( J F/ 2 ,L 2 ,a; 2 ,a,6). 

The proof is analogous to that given in Example [TJ We write the operator — p(+) in 
canonical form Ad 2 R + 69 x + (79^ + -D. Then on an arbitrary basis we have 

(Ad 2 fl + Bd x + C3 R + D)^ n = 0. 

It follows via the usual Wronskian argument that A = B = C = D = 0. 

Similarly we note that the operator E^ - ) = (H_(n+g)H + (n) — H + (n — g)H_(n))/(2n+a+6+l) 
is an even polynomial in (2n + a + b + 1), as is (£ n — r] n )/(2n + a + 6 + 1). Also it is polynomial 
in E 2 and w 2 . Thus H^ - '' = p('(H 2 , L2, w 2 , a, fe) is a symmetry operator which is a polynomial 
function of all of its variables. 
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Now we can compute the structure relations explicitly by evaluating the operators on an 
eigenfunction basis. As we have demonstrated, these relations must then hold everywhere. The 
results are: 

[L 2 ,L 4 ] = -4k 2 qL 3 -4k 2 q 2 L 4 , 

[L 2 , L 3 ] = 2q{L 2 , L 4 } + 4k 2 q 2 L 3 + 8k 2 q 3 L 4 , 

[L 3 ,L 4 ] = 2qL 2 - 2PH (H 2 , L 2 ,u 2 , a, b) , 

6k 2 L\ + {L 2 , L 4 , L 4 } + 6k 2 q{L 3 , L 4 } + 28k 2 q 2 L 2 4 - 4fc 2 gf ( ~ ) {H 2 , L 2 , u 2 ,a, b) 
- 12k 2 pW (H 2 , L 2 ,u 2 , a, b) = 0. 

Here, {A, B} = AB + BA and {A,B,C} is the analogous 6-term symmetrizer of 3 operators. 
A more transparent realization, for R = —4k 2 qL 3 — Ak 2 q 2 L 4 , is 

[L 2 , Li] = R, 

[L 2 ,R] = -8k 2 q 2 {L 2 ,L 4 } - 16k A q A L 4 , 

[Li, R] = 8k 2 q 2 L 2 4 - 8fc 2 gPH (H 2 , L 2 ,u 2 ,a, b) , 

R 2 + 22k 2 q 2 L\ + {L 2 , L 4 , L 4 } - 4fc 2 gP (_) (H 2 , L 2 ,uj 2 , a, b) 



8k 2 q 2 

- 12k 2 P {+) (H 2 , L 2 ,u 2 , a, b) = 0. 
From this realization we see that the symmetries H, L 2 , L 4 generate a closed symmetry algebra. 

Example 6. We consider the TTW system with k = p = q = 1. This is essentially the same 
as Example [5] but with different assumptions. We choose the operators X = d 2 — uj 2 + Ai/x 2 , 
Y = d 2 —uj 2 +A 2 /y 2 and L 2 = M from the previous example where now A\ = \ — a 2 ,A 2 = \ — b 2 . 
Then 

"+/ = Jn K A,mfi = Jn K A,mf 

where A = 2n + a + b + 1, and we form the usual combinations 

L 3 = ~ ++ ~„ = ~1{L 2 ,X- Y } + l-(a 2 - b 2 )H - ±-(X - Y) + ±-[M,X] 

and 

L 1 = i(= + -=_) = i[X,M]. 

From these expressions we deduce that [L 2 ,L 4 ] = — 4(L 3 + L 4 ), as expected. Note that it follows 
from these equations that the second order operator X — Y is a symmetry. Indeed, taking the 
commutator of H with L 3 we find [H, X — Y]L 2 + L 2 [H, X — Y] =0. From this and formal 
adjoint properties one can conclude that [H, X — Y] = 0, so X — Y is a 2nd order symmetry. 
However, Y — Y doesn't belong to the algebra we have already found. In other words our 
standard procedure didn't find the lowest order generator for the symmetry algebra in this case; 
it determined a proper subalgebra of the full symmetry structure. To find the full algebra we 
need to check that the symmetry algebra generated by X — Y, L 2 , H closes at finite order, 
though in this particular example we already know this to be the case. In the next subsection 
we will give another approach to this problem that describes how the missing symmetry can 
always be found and shows that it expressible in terms of the fundamental raising and lowering 
operators. 
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5.1 The symmetry L5 

We investigate the fact that, as shown in Example El our method doesn't always give generators 
of minimal order. For the case k = 1 we know from the structure theory of 2nd order superinte- 
grable 2D systems with nondegenerate potential \28\ [3] that the space of 3rd order symmetries 
is 1-dimensional. Thus, we know that there is a 2nd order symmetry operator L5 for this case, 
independent of L2, H, such that [L2, £5] = L4. We will show how to obtain this symmetry 
from the raising and lowering operators 3±, without making use of multiseparability. Thus, for 
general rational k we look for a symmetry operator L5 such that [£2,-^5] = £4. Applying this 
condition to a formal eigenbasis of functions Vl/ n we obtain the result 

L 2 (L^ n ) = -k\2n + a + b+ l) 2 (L 5 *n) + - ^ - ' "~ - fr n . 

2n + a + + 1 

The general solution is 

u- 1 



4qk 2 \(2n + q + a + b + l)(2n + a + b + 1) 

+ (2n-q + a + b+l)(2n + a + b+l)) + ^ (5 ' 9) 

where f3 n is a rational scalar function. It is easy to check that the quantity in parentheses 
is a rational scalar function of (2n + a + b + l) 2 . Thus we will have a true constant of the 
motion, polynomial in the momenta, provided we can choose f3 n such that the full quantity (15. 9h 
is polynomial in (2n + a + b + l) 2 . To determine the possibilities we need to investigate the 
singularities of this solution at n = — a +H~i ; n = _±9±2Lti?±i > jt j s easy to check that the 
operator has a removable singularity at n = — a+ ^ +1 , In the special case k = 1 the singularities 

at n = — ±1+a + b+l are removable, provided we set j3 n = — ^2(2n+a+b+2)(2n+a+b) ■ Further, the 
expression for L5 in this case k = 1 implies 

{L 5 , L 2 } = -2L 5 + i(L 3 + U) + ^ (a 2 - 6 2 ) . 

Similar calculations show that the symmetries L2, L§, H generate the full closed symmetry 
algebra, which contains our original symmetry algebra properly. For general rational k we can 
set 

a _ Q(H) 

Pn 



4qk 2 (2n + q + a + b + l)(2n - q + a + b + 1) 



and determine a polynomial Q(H) such that the operator has a removable singularity at n = 
— (q + a + b + l)/2, i.e., such that the residue is 0. (Since the solution is a polynomial in 
(2n + a + 6 + l) 2 we don't have to worry about the singularity at n = —{—q + a + b + l)/2.) 
Fixing we can consider the operators as depending on A alone. Note that 



K + - -— - -— 7+ - - (n 2 - h 2 ) K+ - K~ J- - 7+ 



Setting u = —(a + b + l)/2, we see that we have to evaluate the product _ g / 2 _ /x , which 
factors as 



Thus to compute the residue we can pair up terms using the following consequences of the 
recurrence relations above: 

Jn+i J n = J\~ a ~b-2 J n ~ 4(n + 1) (n + a + 6 + 1) (n + a + 1) (n + 6 + 1) , 
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Kl^Kt = K+ A „ K+ ~-^[±L + i + A)[^--l-A 



H H A f H 



A + 2 A ~ 4 v 2^ y V 2^ 

For example, consider the y-dependent terms with p odd. Then the central term in the If -factors 
is multiplication by the constant = —Q. Note that the operators on either side of K^ x pair 
up to give K±K^ 3 which acts as multiplication by a constant. Then we consider the next pair 
KJKZ^, and so on to evaluate the product. If p is even there is no central term and we start 
from the inside by pairing KqK_ 2 , followed by K 2 K^ 4 , and so on. The same procedure works 
for the J- factors, except that the step size is 1, rather than 2. 

The full computation is simple in principle, but technical. It breaks up into 3 cases: 



p even, q odd: 



g m -(.'-*)n((- 1 )(£ + H(£-i 

[ -(-a - b + 2s) (a + b + 2s)(a - b + 2s)(-a + b + 2s) J 

s=l ^ ' 



(?-l)/2 

X 



p odd, g odd: 



(ff-l)/2 . v 

x (-(-a-b + 2s)(a + b + 2s)(a-b + 2s)(-a + b + 2s) J , 



p odd, g even: 



(j>-l)/2 



9/2-1 



( t(-9 - a - 6 + 2/i + l)(-g + a + 6 + 2/i + l)(-g + a - 6 + 2/i + 1) 



X n V 4 

x (-g- a + 6 + 2/i + 1) J . 

Again, by construction one can check that, for general rational k, the operators L2, L$, H 
generate a symmetry algebra that properly contains our original algebra and that closes. The 
new structure relations are somewhat more complicated than before, but computable. Our basic 
point is that the full set of symmetries is generated by the fundamental raising and lowering 
operators E±, even though they themselves are not polynomial symmetries. 



6 One-variable models 

The recurrence operators that we have introduced via special function theory lead almost im- 
mediately to one-variable function space models that represent the symmetry algebra action in 
terms of difference operators on spaces of polynomials. We will illustrate this for two systems, 
the caged anisotropic oscillator and the TTW system, and will describe some of the information 
that can be gleaned from such models. 
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The caged anisotropic oscillator. Prom expressions ()4.5p . (|4.6p and the eigenvalue equa- 
tions for H and L\ it is easy to write down a function space model for irreducible representations 
of the symmetry algebra. Note that since H commutes with all elements of the algebra, it corre- 
sponds to multiplication by a constant E in the model. We let the complex variable t correspond 
to n. Then the action of the symmetry algebra on the space of polynomials f(t) is given by 
difference operators 

$+/(*) = (-4 m )3(4 M2 f (t + l) q {-u + kt — a 2 ) p f(t + q), (6.1) 
$"/(*) = (-4 W )«(4 M2 f (-t - ai) q (-u -kt + l) p f(t - q), (6.2) 
Ef(t) = -2fi(qu + pai + P + qa 2 + q)f(t), L 1 f(t) = -2fi 1 (2t + a 1 + l)f(t), (6.3) 
Mi = PM> M2 = tfM- 
These operators satisfy relations 

$+$- = P 1 (E,L 1 ), *-*+ = fb(£7,Li), 
[L 2 , $ + ] = -4pg/i$ + , [L 2 , $~] = 4pg/x$~. 

(Note that, by using a similarity transformation via the Mellin transform and its inverse, we 
could also transform the model (|6. 1 [) — f)6.3[) into a realization by differential operators in one 
variable.) For a finite dimensional irreducible representation the unnormalized eigenfunctions 
of L 2 are delta functions *&t N (t) = S(t + tjv)- To derive an inner product on the representation 
space we can require that the adjoint of & + is and that L 2 is self-adjoint. 

We can use the model to find several families of finite dimensional and infinite dimensional 
irreducible representations of the symmetry algebra, only some of which correspond to those 
that arise from the quantum mechanical eigenvalue associated with real physical systems. As an 
example of the use of the model to construct finite dimensional representations let us assume 
that ai , a 2 are real and look for a realization such that t at = t° + Nq where N is an integer and 
< t° < q is fixed. Let Nq be the smallest integer such that ^tv = <5(*+*iO is an eigenfunction, 
and let N\ be the largest such integer. Then we must have ^"^'tvo = 0, <I )+x I / 7v 1 = 0. We see 
from the model that one way to accomplish this is to choose 

u - k[t° + N q] + 1 = -po, t° + N lQ + 1 = -g , 

where po, qo are integers such that < po < p and < qo < q. Now set M = N\ — Nq, so that 
the dimension of the representation is M + 1. Simple algebra gives the eigenvalue 

E = 2Mpq/i - 2ji\p(ai - q ) + q(a 2 - Po)] 
for H and eigenvalues 

4/ipg(JVi - N) — 2fi P ( ai -2q -l), N = N , N + 1, . . . , N x , 
for L\. 

The TTW potential. We use expressions (I2.2p . (|2.3p and the eigenvalue equations for H 
and L 2 to write define a function space model for irreducible of the symmetry algebra. The 
results can be presented in a simpler form if we use a gauge transformation <3? n = g(n)*$> n where 

g{n) = (n-q + l) q 

g(n — q) (—u — k(n + a + b + l)) p (—n — a) q 

and introduce a complex variable s, corresponding to n + (a + b + l)/2. Then the action of the 
symmetry algebra on the space of polynomials f(s) is given by difference operators 



L 2 $(s) = -4/e 2 s 2 $(s), E = -2uj (2u + 1 + k(a + b + 1)) , 
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r \ / / a-6 + l\ / a + b + 1 
L 4 $(s) = s + [s + 



From these expressions it is straightforward to show that the space of polynomials in the va- 
riable s 2 is invariant under the action of the one-variable difference operators. For a finite 
dimensional irreducible representation the unnormalized eigenfunctions of L2 are delta functions 
^s N (t) = 8(s + sat). To derive an inner product on the representation space we can require 
that Li and L4 are self- adjoint. 

We can use the difference model to find several families of finite dimensional and infinite 
dimensional irreducible representations of the symmetry algebra, only some of which correspond 
to those that arise from the quantum mechanical eigenvalue associated with real physical sys- 
tems. As an example of the use of the model to construct finite dimensional representations let 
us assume that a, b, 00 are real and look for a realization such that sn = s° + Nq where A is an 
integer and < s° < q is fixed. Let Nq be the smallest integer such that $>Ar = S(s + sn ) is 
an eigenfunction, and let Aq be the largest such integer. We see from the model that one way 
to accomplish this is to choose 

a + b+l\ n „ T a + b + 1 

o = -<?o, 



u + k (s - N q + Q + ^ +1 j + 1 = -po, s° + N t q + 



where po, qo are integers such that < po < P and < qo < q. Now set M = Aq — Nq, so that 
the dimension of the representation is M + 1. We find 

E = 2uj (2pM + a + b + 2kq + 2p + 2) 

for H and eigenvalues 



a + b+1 



2 



-4k 2 ((Aq - N)q + - ' 2 ' - + g j , A = N , A + 1, . . . , Aq, 
for Li- 

The system (13.11) . For completeness we give the simple one variable model for this case: 

L 2 F(n) = n 2 F{n), $+F(n) = p p (F{n + P ) + (-i) p+cl F(n - P )) , 

= (F(fi +p) - {-If^F^l - p)) . 



7 Stackel transforms and the recurrence method 

The theory of Stackel transforms does not guarantee that if two classical 2D superintegrable 
systems are related by a Stackel transform and if one system is quantum superintegrable, then 
the other system is also quantum superintegrable and the quantum systems are related by 
a Stackel transform. Some additional conditions must be fulfilled, [18]. However, if one of the 
systems is known to be quantum superintegrable via our recurrence relation method, then it 
is automatic that the second quantum system is also superintegrable and a quantum Stackel 
transform of the first. We will give a single example which makes clear the general proof. 

In |16j there was introduced a new family of Hamiltonians with a deformed Kepler-Coulomb 
potential dependent on an indexing parameter k which was shown to be related to the TTW 
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oscillator system system via coupling constant metamorphosis. The authors showed that this 
deformed Kepler system is classically superintegrable for all rational k = p/q, and in [22J we 
used the canonical equations for higher order symmetry operators to show that that it is quan- 
tum superintegrable. Here we use the recurrence relations obeyed by the eigenfunctions of the 
separating symmetry operators to give a new proof of quantum superintegrability and, also, to 
obtain the structure equations. 

As stated above, expressed in polar coordinates r, 6, the quantum TTW system is H^> = 

or 

d r + -Or + \d$ - u,V + - a + 2 { ) * = E9. (7.1) 
r t 1 r 2 sm(k9) r jS cos z (k6)J 

The deformed Kepler-Coulomb system, expressed in polar coordinates R, <f>, is H'^f = £^ or 

8 « + it 8 " + w 8 i ~ I + m J ikm + w Jwm )* = £ *- (7 - 2) 

Now note that if we divide both sides of expression (|7.1|) by r 2 , rearrange terms, and make 
the change of variables R = r 2 , <f> = 26, then this expression is identical to (|7.2p with the 
identifications 

I o Z. (7.3) 

Thus (|7.2p is a Stackel transform of (|7.ip and the two systems are Stackel equivalent, [TBI [T6[ [22] . 

The principal observation that we need to make is that both systems have exactly the same 
formal eigenfunctions (|5.ip . modulo variable substitution and identifications (17.3p . and exactly 
the same separation equations and recurrence formulas. Thus, substituting into (|5.2p we see 
that the energy levels for the deformed Kepler-Coulomb system are 

£ 



(2(m + nk) + 1 + (a + b + l)k) 



2 ' 



In the expressions for the raising and lowering operators of the TTW system we replace the 
operator E ~ H by the constant AZ to get the the raising and lowering operators for the 
deformed Kepler-Coulomb system. Similarly the structure equations for the Kepler-Coulomb 
system are obtained by simple permutations 

u? O AH', H o AZ. 

We find 

[L2, L4] = R, 

[L 2 ,R] = -8k 2 q 2 {L 2 ,L 4 } - 16k A q A L 4 , 

[Li, R] = Sk 2 q 2 L\ - 8k 2 qP^ (WZ 2 , L 2 , AH', a, b) , 

^pR 2 + 22k 2 q 2 L\ + {L 2 , L 4 , L 4 } - Ak 2 qP { ~ ] (16Z 2 , L 2 , AH', a, b) 
- l2k 2 P i+) (WZ 2 , L 2 ,AH', a, b) = 0. 



Similarly the operator L5 can be added to the system. 
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8 Discussion and conclusions 

We have developed a new method for verifying quantum superintegrability for 2D systems and 
given several applications to families of such systems, most notably the caged anisotropic os- 
cillator, the Tremblay-Turbiner-Winternitz system and the deformed Kepler-Coulomb system 
and given new proofs of superintegrability for all rational k. The method relies on the assump- 
tion that the system has a second order symmetry operator, so that the Schrodinger eigenvalue 
equation H^> = separates in a set of orthogonal coordinates u\, U2 determined by the 
symmetry and so that the separated eigenfunctions ^ = Ui{u{)U2{u2) satisfy computable recur- 
rence relations. In practice this means that the separated solutions need to be of hyper geometric 
type. Then one employs the recurrence relations to construct operators that commute with H 
on a formal eigenbasis. We used our earlier developed canonical form for higher (> 2) order 
symmetry operators to show that operators commuting with H on a formal eigenbasis must be 
actual symmetry operators, further that operator identities verified on formal eigenbases must 
hold identically. Using this approach one can obtain explicit, though complicated, expressions 
for the higher order symmetry operators. 

We saw that in the case of the TTW potential our method didn't lead immediately to the 
lowest order generators, but that they could be found and expressed in terms of our raising and 
lowering operators. We have no proof as yet that we have found the lowest order generator for 
all rational k but this is the case for all examples that we know. Provable determination of 
the maximal symmetry algebra is a topic for future research. By acting on formal eigenbases 
we were able to compute symmetry algebras for each of the systems we studied and show that 
these algebras were closed under commutation. One striking result of these computations was 
that the structure equations for the symmetries were quite explicit and much simpler that the 
expressions for the symmetry operators themselves. In essence, one can determine the structure 
equations without knowing the higher order generating symmetries! 

In each case the action of the symmetries on a formal eigenbasis led us to a simple model 
of the associated symmetry algebra and its representations, in terms of difference operators 
in one variable. This greatly simplifies the analysis of the structure equations, classification of 
irreducible representations of the symmetry algebra and determination of the spectral properties 
of the generating symmetries. Of course adding the generator L5 to the algebra complicates these 
models, though by construction it can always be realized as a difference operator. This is an 
issue for future research. Such models have independent interest [Tj IH1 125]. 

Finally, we gave an example of the use of the Stackel transform to map one superintegrable 
family of systems to another, while preserving the structure equations. In this case we mapped 
the TTW system to a deformed Kepler-Coulomb system and determined, for the first time, the 
structure equation for the Kepler-Coulomb system. 

It is interesting that this approach to quantum superintegrability with structure results for 
symmetry algebras has preceded the classical approach; usually the reverse is true. To our 
knowledge, the closure of classical symmetry algebras has not been proven for the systems 
considered here. We are actively investigating the classical analog of the quantum construction. 
It is clear that these methods have much greater applicability than just the examples treated 
here. Indeed, all of the systems studied in |22] could be so analyzed. Further all our methods 
can clearly be extended to systems in dimensions > 2. There appears to be no obstacle other 
than growing complexity. 

When it applies, the recurrence relation method is much simpler than our earlier introduced 
canonical form method for verification of quantum superintegrability and it provides much 
more information, including the structure equations. However, the recurrence method requires 
a detailed knowledge of recurrence relations obeyed by the separated solutions and there are 
examples of superintegrable systems where no such relations appear to exist [301 Section 4.2]. 
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Thus the canonical form method for higher order symmetry operators seems to be more general. 
Furthermore, we used the canonical form at a crucial point in the recurrence approach to show 
that computations valid on formal eigenbases actually held identically. The approaches are 
mutually complementary. 
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